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A study of the satellite 1962/3a (Alouette) during the period April 27, 
1963 to March 8, 1964 revealed a long-period perturbation (approximately 
140 days, the period of the argument of perigee) in the eccentricity of its 
orbit, plus shorter -period variations. An analysis of the possible sources 
of these perturbations was undertaken, neglecting air drag since perigee 
was about 1000 km. 

Mean orbital elements for this period (Table 1), computed using 
Brouwer's satellite theory (1), were obtained from the Goddard computing 
center. These elements take into account the effects due to zonal har- 
monics (through J s ) in the earth's gravitational potential. Figure 1 
illustrates the variation of e with time. fyu ’i 


LUNAR AND SOLAR GRAVITATIONAL EFFECTS 

The disturbing function for both the sun and moon may be written in 
the form 
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where 


G = gravitational constant 

r = geocentric distance of satellite 

r' = geocentric distance of sun (or moon) 

S 2 , S 3 ,... = Legendre polynomials 

# _ mass of disturbing body 

mass of disturbing body + mass of earth 

The variation equation for eccentricity is 

de _ V 1-e 2 BR 1-e 2 BR 

na 2 e ^ na 2 e 

where 

n = mean motion of satellite 
a = semi-major axis of satellite’s orbit 
to = argument of perigee 
M = mean anomaly of satellite. 

Assuming mean values for the elements and considering only the long- 
period terms from the second-order part of the disturbing function 
[Kozai (2)] , a first approximation to Se reveals no terms with ampli- 
tudes greater than 0(10“ 7 ). Hence, luni-solar gravitational effects were 
considered negligible. 

SOLAR RADIATION PRESSURE 

Neglecting the effect of the earth’s shadow, the solar radiation pres- 
sure disturbing function is [Kaula (3)] 
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-—sin i sin e cos (a>+A. 0 ) + — sin i sin e cos (oj-\ 0 ) 

2 2 


where 



i = inclination of satellite's orbit plane to earth's equatorial plane 
H = longitude of ascending node 
e = obliquity of the ecliptic 
= mean longitude of the sun 
F = force constant = [-4.63X10 -5 dynes/cm 2 ] X A 
A = presentation area of the satellite* 1.86 meters 2 
m = mass of satellite = 145.15 kg. 

(As an approximation, the satellite's "presentation area" was assumed 
to be 1/4 its surface area). 


Substituting into the variation equation and again using mean values 
for the elements, a first approximation yields 


Se = .000040 cos (oj-Q+Xq) + .000007 cos (co+Q-K 0 ) 

- .000007 cos (<y+A. 0 ) + .000003 cos (co-\ Q ). 

Values of e corrected for solar radiation pressure (labeled e c ) are given 
in Table 1. Figure 2 gives a plot of these values vs. time. 
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LEAST SQUARES FIT AND HARMONIC ANALYSIS 


A least squares fit of the type 

e = e Q + A sin co + B cos 2 co 

was applied to the data in Figure 2 - the resulting equation was 

e = .002482 + .000161 sin co + .000133 cos 2 co. 

Figure 3 indicates the closeness of the fit — residuals are given in 
Table 1. The standard deviation was .000036. 

In addition, a 12-point harmonic analysis was performed on the data 
in Figure 2 for one period of the argument of perigee, beginning 99 days 
after April 27. The series obtained was 

e = .002490 + .000149 sin + .000117 cos 2 co — ■ 000027 sin 3 co 

- .000027 cos 4 co - .000011 sin 5 co. 


DISCUSSION 

It was found that the mean values of the eccentricity of Alouette's 
orbit listed in Table 1 can be represented by the series 

e = .002482 + .000161 sin a> + .000133 cos 2 co + .000040 cos ( co — 0 + 

+ .000007 cos (ai+fi-X.©) - .000007 cos (ca+X.©) + .000003 cos (w-X.©). 

A more accurate representation for a 140 day period beginning 99 days 
after April 27 is 
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e = .002490 + .000149 sin a> + .000117 cos 2 w - .000027 sin 3 oj - .000027 cos 4 a; 


- .000011 sin 5 co + .000040 cos (aj-ft+^ 0 ) + .000007 cos (a>+ft-^ 0 ) 

- .000007 cos (co+k Q ) + .000003 cos (oo~k Q ). 

• In addition, an analysis was made of the argument of perigee. Solar re- 

sonances were found for terms in the solar gravitational disturbing function 
with arguments 2 ( k Q - 2ft© + ft ) , k Q - o> 0 - 2 ft 0 + ft , and k Q + co Q - 2ft 0 + ft , 
so the effects of these terms could not be computed by the method used 
previously. However, corrections were made for luni- solar secular 
effects and solar radiation pressure. Precession and nutation effects 
were computed using results obtained by Kozai (4), and were found to be 
relatively insignificant. Then the secular part of co was removed by 
means of a least squares fit of the type 

= co 0 + CO (t-t 0 ), 

and the residuals were plotted against time (Figure 4). It is interesting 
to note that the remaining variation appears similar to that found in the 
eccentricity (Figure 2), the period again approximating that of the argu- 
ment of perigee. 
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Table 1 

Mean Values of co and e — e Corrected for Solar Radiation 
Pressure and Comparison with Least Squares Fit 


1 - to 

0) 

e 

e 

c 

Residuals 

0 

146.165 

.00262 

' ' ' 1 

+2 

7 

129.203 

.00259 


+ 3 

21 

92.246 

.00248 

.00250 

-1 

28 

73.179 

.00249 

.00251 

-2 

43 

35.239 

.00256 

.00260 

-2 

50 

19.141 

.00253 

.00257 

-7 

57 

2.565 

.00251 

.00256 

-6 

64 

346.096 

.00247 

.00251 

-5 

71 

328.431 

.00239 

.00244 

-2 

78 

311.372 

.00230 

.00235 

+ 1 

85 

289.665 

.00210 

.00214 

-9 

99 

246.193 

.00220 

.00224 

-1 

106 

225.570 

.00237 

.00241 

+ 5 

113 

207.545 

.00244 

.00248 

0 

127 

174.010 

.00256 

.00260 

-3 

134 

159.097 

.00258 

.00262 

-2 

148 

126.361 

.00255 

.00259 

+2 

162 

89.885 

.00248 

.00250 

-1 

169 

67.343 

.00253 

.00255 

+ 1 

176 

51.490 

.00254 

.00256 

-2 

183 

34.968 

.00259 

.00260 

-2 

190 

18.434 

.00261 

.00262 

-2 

204 

346.682 

.00255 

.00256 

0 

211 

329.351 

.00244 

.00245 

-1 

218 

311.139 

.00237 

.00237 

+ 3 

225 

292.302 

.00223 

.00222 

-2 

232 

267.748 

.00218 

.00217 

-2 

239 

246.464 

.00230 

.00228 

+4 

246 

227.052 

.00243 

.00241 

+6 

253 

209.096 

.00255 

.00253 

+6 

260 

191.239 

.00264 

.00261 

+4 

267 

176.018 

.00268 

.00265 

+ 3 

274 

160.457 

.00270 

.00268 

+4 

281 

145.475 

.00273 

.00271 

+9 

288 

127.954 

.00261 

.00259 

+ 1 

295 

110.000 

.00253 

.00251 

-2 

302 

90.716 

.00253 

.00251 

0 

309 

70.641 

.00256 

.00254 

+ 1 

316 

53.231 

.00260 

.00257 

0 


a = 1.1589 earth radii i = 80.466° 


t - t 0 is the number of days since April 27, 1963. 


os is measured in degrees. 
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Figure 1. Eccentricity of 1962 /3a (Published values) 
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Figure 2. Eccentricty of 1962 /3a (Published values less solar radiation pressure) 
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Figure 3. Eccentricity of 1962 /3a 
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